Calculating the axes of rotation for the subtalar and talocrural joints using 3D bone reconstructions

Orientation of the subtalar joint axis dictates inversion and eversion movements of the foot and has been the focus of evolutionary and clinical studies for a number of years. Previous studies have measured the subtalar joint axis against the axis of the whole foot, the talocrural joint axis and, recently, the principal axes of the talus. The present study introduces a new method for estimating average joint axes from 3D reconstructions of bones and applies the method to the talus to calculate the subtalar and talocrural joint axes. The study also assesses the validity of the principal axes as a reference coordinate system against which to measure the subtalar joint axis. In order to define the angle of the subtalar joint axis relative to that of another axis in the talus, we suggest measuring the subtalar joint axis against the talocrural joint axis. We present corresponding 3D vector angles calculated from a modern human skeletal sample. This method is applicable to virtual 3D models acquired through surface-scanning of disarticulated ‘dry’ osteological samples, as well as to 3D models created from CT or MRI scans

Human Remains from the Pleistocene-Holocene Transition of Southwest China Suggest a Complex Evolutionary History for East Asians

BACKGROUND: Later Pleistocene human evolution in East Asia remains poorly understood owing to a scarcity of well described, reliably classified and accurately dated fossils. Southwest China has been identified from genetic research as a hotspot of human diversity, containing ancient mtDNA and Y-DNA lineages, and has yielded a number of human remains thought to derive from Pleistocene deposits. We have prepared, reconstructed, described and dated a new partial skull from a consolidated sediment block collected in 1979 from the site of Longlin Cave (Guangxi Province). We also undertook new excavations at Maludong (Yunnan Province) to clarify the stratigraphy and dating of a large sample of mostly undescribed human remains from the site. METHODOLOGY/PRINCIPAL FINDINGS: We undertook a detailed comparison of cranial, including a virtual endocast for the Maludong calotte, mandibular and dental remains from these two localities. Both samples probably derive from the same population, exhibiting an unusual mixture of modern human traits, characters probably plesiomorphic for later Homo, and some unusual features. We dated charcoal with AMS radiocarbon dating and speleothem with the Uranium-series technique and the results show both samples to be from the Pleistocene-Holocene transition: ∼14.3-11.5 ka. CONCLUSIONS/SIGNIFICANCE: Our analysis suggests two plausible explanations for the morphology sampled at Longlin Cave and Maludong. First, it may represent a late-surviving archaic population, perhaps paralleling the situation seen in North Africa as indicated by remains from Dar-es-Soltane and Temara, and maybe also in southern China at Zhirendong. Alternatively, East Asia may have been colonised during multiple waves during the Pleistocene, with the Longlin-Maludong morphology possibly reflecting deep population substructure in Africa prior to modern humans dispersing into Eurasia

Supraorbital morphology and social dynamics in human evolution

Uniquely, with respect to Middle Pleistocene hominins, anatomically modern humans do not possess marked browridges, and have a more vertical forehead with mobile eyebrows that play a key role in social signalling and communication. The presence and variability of browridges in archaic Homo and their absence in ourselves have led to debate concerning their morphogenesis and function, with two main hypotheses being put forward; that browridge morphology is the result of the spatial relationship between the orbits and the braincase, and that browridge morphology is significantly impacted by biting mechanics. Here we virtually manipulate browridge morphology of an archaic hominin (Kabwe 1), showing that it is much larger than the minimum required to fulfil spatial demands and that browridge size has little impact on mechanical performance during biting. Since browridge morphology in this fossil is not driven by spatial and mechanical requirements alone, the role of the supraorbital region in social communication is a potentially significant factor. We propose that conversion of the large browridges of our immediate ancestors to a more vertical frontal in modern humans allowed highly mobile eyebrows to display subtle affiliative emotions

Alpha shapes: Determining 3D shape complexity across morphologically diverse structures

Background. Following recent advances in bioimaging, high-resolution 3D models of biological structures are now generated rapidly and at low-cost. To utilise this data to address evolutionary and ecological questions, an array of tools has been developed to conduct 3D shape analysis and quantify topographic complexity. Here we focus particularly on shape techniques applied to irregular-shaped objects lacking clear homologous landmarks, and propose the new ‘alpha-shapes’ method for quantifying 3D shape complexity. Methods. We apply alpha-shapes to quantify shape complexity in the mammalian baculum as an example of a morphologically disparate structure. Micro- computed-tomography (μCT) scans of bacula were conducted. Bacula were binarised and converted into point clouds. Following application of a scaling factor to account for absolute differences in size, a suite of alpha-shapes was fitted to each specimen. An alpha shape is a formed from a subcomplex of the Delaunay triangulation of a given set of points, and ranges in refinement from a very coarse mesh (approximating convex hulls) to a very fine fit. ‘Optimal’ alpha was defined as the degree of refinement necessary in order for alpha-shape volume to equal CT voxel volume, and was taken as a metric of overall shape ‘complexity’. Results Our results show that alpha-shapes can be used to quantify interspecific variation in shape ‘complexity’ within biological structures of disparate geometry. The ‘stepped’ nature of alpha curves is informative with regards to the contribution of specific morphological features to overall shape ‘complexity’. Alpha-shapes agrees with other measures of topographic complexity (dissection index, Dirichlet normal energy) in identifying ursid bacula as having low shape complexity. However, alpha-shapes estimates mustelid bacula as possessing the highest topographic complexity, contrasting with other shape metrics. 3D fractal dimension is found to be an inappropriate metric of complexity when applied to bacula. Conclusions. The alpha-shapes methodology can be used to calculate ‘optimal’ alpha refinement as a proxy for shape ‘complexity’ without identifying landmarks. The implementation of alpha-shapes is straightforward, and is automated to process large datasets quickly. Beyond genital shape, we consider the alpha-shapes technique to hold considerable promise for new applications across evolutionary, ecological and palaeoecological disciplines

3D bite modeling and feeding mechanics of the largest living amphibian, the Chinese Giant Salamander Andrias davidianus (Amphibia:Urodela)

Biting is an integral feature of the feeding mechanism for aquatic and terrestrial salamanders to capture, fix or immobilize elusive or struggling prey. However, little information is available on how it works and the functional implications of this biting system in amphibians although such approaches might be essential to understand feeding systems performed by early tetrapods. Herein, the skull biomechanics of the Chinese giant salamander, Andrias davidianus is investigated using 3D finite element analysis. The results reveal that the prey contact position is crucial for the structural performance of the skull, which is probably related to the lack of a bony bridge between the posterior end of the maxilla and the anterior quadrato-squamosal region. Giant salamanders perform asymmetrical strikes. These strikes are unusual and specialized behavior but might indeed be beneficial in such sit-and-wait or ambush-predators to capture laterally approaching prey. However, once captured by an asymmetrical strike, large, elusive and struggling prey have to be brought to the anterior jaw region to be subdued by a strong bite. Given their basal position within extant salamanders and theirPeer ReviewedPostprint (published version

Finite element analysis of a micromechanical model of bone and a new 3D approach to validation

Finite Element Analysis (FEA) is now widely used to analyse the mechanical behaviour of bone structures. Ideally, simulations are validated against experimental data. To date, validation of Finite Element Models (FEMs) has been 2 Dimensional (2D) only, being based on comparison with surface-mounted strain gauge readings. In this study we present a novel 3-Dimensional (3D) approach to validation that allows comparison of modelled with experimental results between any two points in 3D space throughout the structure, providing magnitude and direction data for comparison, internally and externally. Specifically, we validate a FEM of a rat tibia, including trabecular network geometry, using a material testing stage housed within a microCT scanner. We further apply novel landmark based morphometric approaches to more effectively compare modelled and experimental results. 542 landmark points on the cortical and trabecular bone surfaces of the model were selected and validated in 3D against experimental data. This approach may hold considerable potential in fields wherein a better understanding of the mechanical behaviour of trabecular networks is important, e.g., the studies of osteoporosis and trabecular loss after orthopaedic implant insertion. © 2012 Elsevier Ltd

Finite element micro-modelling of a human ankle bone reveals the importance of the trabecular network to mechanical performance: New methods for the generation and comparison of 3D models

Most modelling of whole bones does not incorporate trabecular geometry and treats bone as a solid non-porous structure. Some studies have modelled trabecular networks in isolation. One study has modelled the performance of whole human bones incorporating trabeculae, although this required considerable computer resources and purpose-written code. The difference between mechanical behaviour in models that incorporate trabecular geometry and non-porous models has not been explored. The ability to easily model trabecular networks may shed light on the mechanical consequences of bone loss in osteoporosis and remodelling after implant insertion. Here we present a Finite Element Analysis (FEA) of a human ankle bone that includes trabecular network geometry. We compare results from this model with results from non-porous models and introduce protocols achievable on desktop computers using widely available softwares. Our findings show that models including trabecular geometry are considerably stiffer than non-porous whole bone models wherein the non-cortical component has the same mass as the trabecular network, suggesting inclusion of trabecular geometry is desirable. We further present new methods for the construction and analysis of 3D models permitting: (1) construction of multi-property, non-porous models wherein cortical layer thickness can be manipulated; (2) maintenance of the same triangle network for the outer cortical bone surface in both 3D reconstruction and non-porous models allowing exact replication of load and restraint cases; and (3) creation of an internal landmark point grid allowing direct comparison between 3D FE Models (FEMs). © 2012 Elsevier Ltd

Finite element analysis of a micromechanical model of bone and a new 3D approach to validation

Finite Element Analysis (FEA) is now widely used to analyse the mechanical behaviour of bone structures. Ideally, simulations are validated against experimental data. To date, validation of Finite Element Models (FEMs) has been 2 Dimensional (2D) only, being based on comparison with surface-mounted strain gauge readings. In this study we present a novel 3-Dimensional (3D) approach to validation that allows comparison of modelled with experimental results between any two points in 3D space throughout the structure, providing magnitude and direction data for comparison, internally and externally. Specifically, we validate a FEM of a rat tibia, including trabecular network geometry, using a material testing stage housed within a microCT scanner. We further apply novel landmark based morphometric approaches to more effectively compare modelled and experimental results. 542 landmark points on the cortical and trabecular bone surfaces of the model were selected and validated in 3D against experimental data. This approach may hold considerable potential in fields wherein a better understanding of the mechanical behaviour of trabecular networks is important, e.g., the studies of osteoporosis and trabecular loss after orthopaedic implant insertion